arXiv:1506.03952vl [cond-mat.soft] 12Jun2015 


Condensed Matter Physics, 2015, Vol. 18, No 2, 23001: 1 -[i^ 
DOI: 10.5488/CMP.18.23001 
http://www.icmp.lviv.ua/journal 




Anomalous Brownian motion of colloidal particle in a 
nematic environment: effect of the director 
fluctuations 

T. Turi\i®, A. Brodirl®, V.G. NazarenkcP 


^ Institute of Physics of the Nationai Academy of Sciences of Ukraine, 46 Nauky Ave., 03028 Kyiv, Ukraine 
^ Liquid Crystai Institute, Kent State University, 1425 University Espianade, 44242 Kent, OH, USA 
^ Nationai Technicai University of Ukraine "Kyiv Polytechnic Institute", 37 Peremogy Ave., 03056 Kyiv, Ukraine 

Received October 1,2014, in final form December 3, 2014 

As recently reported [Turiv T. et al.. Science, 2013, 342, 1351], fluctuations in the orientation of the liquid 
crystal (LC) director can transfer momentum from the LC to a colloid, such that the diffusion of the colloid 
becomes anomalous on a short time scale. Using video microscopy and single particle tracking, we investigate 
random thermal motion of colloidal particles in a nematic liquid crystal for the time scales shorter than the 
expected time of director fluctuations. At long times, compared to the characteristic time of the nematic director 
relaxation we observe typical anisotropic Brownian motion with the mean square displacement (MSD) linear in 
time T and inversly proportional to the effective viscosity of the nematic medium. At shorter times, however, the 
dynamics is markedly nonlinear with MSD growing more slowly (subdiffusion) or faster (superdiffusion) than t. 
These results are discussed in the context of coupling of colloidal particle's dynamics to the director fluctuation 
dynamics. 
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1. Introduction 

Unlike the visible world around us, the atomic, molecular and small particle worlds are In a state of 
constant motion. This motion is widely occurring in nature and plays important role In physics, chemistry, 
biology and engineering (1). The physical approach to this motion (Brownian motion (2)), developed early 
in the 20th century by Einstein (3l|4), Smoluchowski O, and Langevin (6), stih forms the basis for our 
understanding of these stochastic dynamics, with the main result, derived by Einstein (SHU, being that 
the mean squared displacement (MSD) (Ar^) of a particle undergoing Brownian motion in a Newtonian 
fluid increases linearly with time, (Ar^(T)) = 6 Dt, where D is the diffusion constant, for a spherical 
particle with hydrodynamic radius if In a fluid with viscosity rj, the diffusion coefficient is given by the 
Stokes-Einstein relation D- where is the Boltzmann constant, T is the temperature and C, is 

the viscous friction coefficient, which under no-slip conditions is given by = &nRrj. The normal, linear 
diffusion regime for a particle of mass m Is established on time scales long compared to the microscopic 
time m/(, which is t)q)icaUy in the nanosecond range. These results are valid for Brownian motion under 
the influence of two forces, the viscous frictional force linear in particle velocity, F = -(v, and a random 
force with white-noise-spectrum due to random collisions with surrounding particles. 

Brownian particles in complex systems may, however, exhibit quite different dynamics, reflecting the 
properties of particles themselves (e.g., the effects of a crossover from short-time anisotropic to long-time 
isotropic diffusion for ellipsoidal particles In water (T)) as weU as the local properties of the host medium 
that may be inhomogeneous, exhibit nonlinear friction or elastic properties, etc., with (Ar^(T)> oc t“, 
where the exponent a may be different from 1. For example, colloidal particles in polymer networks 
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(HO, in F-actin networks (TOl, in surfactant formed lyotropic liquid crystal (TT) may exhibit a subdiffusive 
behavior with 0 < a < 1; superdiffusion, a > 1, was observed in concentrated suspensions of swimming 
bacteria (12] [T^ and in “living polymers” (14). 

The behaviour of colloidal particles in a nematic liquid crystal (NLC) is in many respects more com¬ 
plicated as compared to the isotropic fluid host. First of all, the particle sets a certain director distortion 
around itself, due to the anisotropic nature of surface interactions (surface anchoring) (T^. These di¬ 
rector distortions lead to long-range elastic interactions of the particle with the bounding walls. Second, 
the orientational nematic order leads to an anisotropy of the Stokes drag 11611251 . As a result. Brownian 
motion in a nematic host becomes anisotropic with two different diffusion constants and D^, corre¬ 
sponding to the directions parallel and perpendicular to the nematic director n. Third, although the LC 
medium is homogeneous, the average axes of orientation fluctuate in time and space and thus influence 
the diffusive regimes (26|[27j. In addition to normal diffusion, the particle experiences two anomalous 
regimes, with MSB growing slower (subdiffusion) and faster (superdiffusion) than t. The anomalous dif¬ 
fusion occurs at time scales that correspond to the relaxation times of director fluctuations (26). All three 
regimes of diffusion are anisotropic, with the MSB being larger for the motion along the director. Once 
the nematic is melted, the diffusion becomes normal and isotropic. 

In this paper, we elaborate in detail the Brownian dynamics of colloidal particles in a nematic host 
on time scales shorter than the expected time of director fluctuations Tj-eiax- Following the experiments 
described in OS)> we first give a detailed description of the particle tracking experiments with an ex¬ 
tensive analysis of the possible errors in the observed trajectories of the colloids. This is followed by an 
analysis of the diffusion of colloids, which turns out to be nonlinear with time for a certain time lag, as 
it is described in (26). In the theoretical section, we present a model consideration of the diffusion of the 
colloidal particle in LC, where we explain the anomalous regime by coupling of the particle motion and 
director fluctuations. 


2. Experiment 

2.1. Material 

The refractive indices for ordinary (Wq) and extraordinary (Tie) modes of light propagation in a ne¬ 
matic are different from each other. The director n is also the local optic axis. Whenever the director 
fluctuates, birefringence An - Ug- Uo and the associated “lens” effect of the distorted optic axis around 
the particle translate these fluctuations into phantom drifts of the image. To overcome the problem, we 
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Figure 1. (Color online) Wavelength dependence of the refractive indices and optical anisotropy of the 
NLC IS-8200 which chemical structure formula is placed in the plot. 

use the nematic IS-8200 synthesized at Merck (Germany) with ultra-low birefringence An - 0.0015 (at 
light wavelength A = 520 nm) (28). IS-8200 is a thermotropic (solvent free, single-component) material 
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with the nematic state in the temperature range 47-53°C, above which it melts into the isotropic fluid, 
figurej^ Small An suppresses the image shifts caused by fluctuations. 


2.2. Director alignment 

Alignment of the director at both the particle’s surface and the bounding plates influences the diffu¬ 
sion and thus needs to be controlled. We explored the alignment at the surface of spheres to be perpendic¬ 
ular to the surface. In his case, the particles were functionalized with dimethyl-octadecyl- [3(trimethoxysi- 
lyDpropyl] ammonium chloride (DMOAP) (29). The overall uniform orientation of the nematic was set by 
two glass plates covered with rubbed polyimide PI-2555 (Nissan Chemicals) ahgnment layers that pro¬ 
duce a uniaxial planar alignment no = (1,0,0) = const in the cell. The locally distorted director around 
the spheres should smoothly transform into no (22). The resulting equilibrium director configuration is 
of a dipolar type (with a point defect — hyperbolic hedgehog accompanying the sphere, figure]^ for the 
normal anchoring. The director distortions around the particle (30) lead to a repulsion from the bounding 
substrates (M). The particles levitate (^ in the bulk at some height determined by the balance of gravity 
and elastic forces. 




(a) (b) 

Figure 2. (Color online) (a) Schematic director field around a sphere with normal surface anchoring and 
(b) optical image of 5 pm silica sphere of the dipole structure immersed in IS-8200 placed between crossed 
polarizers. 





f(s) 


(b) 


Figure 3. (Color online) (a) Intensity of the transmitted light as a function of the appUed voltage, (b) In¬ 
tensity of the transmitted light as a function of time after voltage is switched off. IS-8200 liquid crystal, 
homeotropic alignment, crossed polarizers. CeU thickness 20 pm. 
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2.3. Time scale 

The director field n (r, t) fluctuates in time and space (32). The relaxation time of fluctuations is much 
slower than the relaxation time of vorticity ~ pl^lrj ~ 1 ps, which is the time needed by the perturbed 
fluid of density p and viscosity p ~ 0.1 Pa ■ s to flow over the distance 1. For a director perturbation with 
a wavevector q, this time is r^eiax ~ where peff and K are the effective viscosity and elastic 

constant. A uniaxial nematic features five independenf viscosities and at least four Frank elastic con¬ 
stants, thus peff and K are complex combinations of these that depend on the director field (33). In a flat 
cell of a finite thickness h, the fluctuation spectrum is restricted; in case of strong director anchoring at 
the boundaries, the minimum wavevector component in the z-direction perpendicular to the cell bound¬ 
aries is -nih (34). We are interested in wavevectors not much larger than q - n/d, where d is the 
sphere’s diameter, since perturbations with wavelength much shorter than d will produce on average 
a negligible effect. The corresponding range of viscous relaxation times is thus < Treiax < t/,, where 
T/i = qeffh^l^^K and = peffd^ To estimate these, we experimentally determined the quantity 
Peff/K in the splay Frederiks transition of lS-8200 and found PefflK =; 10^^ s/m^, figurej^ 

2.4. Optical particle tracking: experimental and data analysis procedures 

Experimental setup consists of a CMOS high-speed video camera MotionBlitz EoSens minil (Microtron 
GmbH) mounted on an inverted microscope Nikon TE2000-U with a 100 x 1.3 N.A. immersion objective. 
The camera is capable of a 5000 fps maximum frame rate (time resolution 0.2 ms); the images are gray¬ 
scale with 256 gradations (bit depth 8). The pixel size is 14 x 14 microns, so that the magnification is 
140 nm/pixel, with the Airy disk covering ~ 100 pixels. Eor temperature control, we used a Linkam LTS120 
heating stage (accuracy 0.3“C). 

for a point light source, its image is diffraction-blurred into an Airy disk, whose intensity profile in the 
image plane has a root-mean-square (r.m.s.) size of s = 0.44AM/2 N.A., where A is the wavelength, N.A. 
the numerical aperture, and M the transverse magnification of the optical system. N.A. should be max¬ 
imized in order to decrease the diffraction blur. Meanwhile, since the image is pixelised as it is formed 
on a detector matrix, larger magnification M implies more pixels in the image, and thus more informa¬ 
tion and better ultimate position accuracy. In practice, N.A. and M are determined by the objective lens, 
typically with N.A. ~ 1 and M ~ 100, so that s ~ 10 pm. 

Particle trajectories are almost invariably analyzed in terms of mean square displacement (33, MSB, 
i.e., (Ax^(t)) and (Ay^fr)), where t is the time lag and angular brackets stand for ensemble average. In 
an isotropic medium (Ajc^(t)) and (Ay^fr)) are equal, whereas in an anisotropic medium, such as NEC, 
they are different. Ensemble average cannot be obtained from a single trace. Instead, one can do a time 
average which is believed to be the same as the ensemble average in the limit of infinite averaging time. 
Specifically, one computes 

{Ax^{T))^{[x[t + T)-x[t)f), (Ay^(T)) = ([y(f-i-T)-y(r)]^), (1) 

where angular brackets now stand for time average. Even though time averages equal the ensemble 
averages in an ergodic system, time averages over a trace of finite length and time step possess specific 
statistical errors (36) . 

The measurement errors {Sxi, Syi) (assumed to have zero mean) in particle coordinates (x,-, y/) yield a 
positive additive contribution to MSB computed through equation Q, as it is a quadratic form. Assuming 
(Xi,yi) = (Xi-i-5xi,yj-i-5y,), where (Xj,y/) are the true coordinates, and denoting x/ = Xf, Sxi - 5xf, where 
t - iAt, one gets for (Ax^(t)) (and similarly for (Ay^(T))) 

(Ax^(t)) = {[[Xt+r + Sxt+r)-{Xt + Sxt)f) 

= {{Xt+r-Xtf) + {[X,+r-X,)[5Xt+r-SXt)) + {[5Xt+T-SXtf), ( 2 ) 

where the first term is the “true” MSB. Assuming that the errors are uncorrelated with the coordinates, 
(xSx) - 0, the second term in equation (3 averages to zero. Assuming further that the errors of different 
trace points are uncorrelated, {Sxt+rSxt) - 0, and have the same variance (Sx^) - Ag, the last term in 
equation (3 yields 2Ag. Thus, 

(Ax2(T)) = (Ax^(T))-t2Ao, (3) 
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Figure 4. (Color online) (a) Grayscale image of 5 ^im particle glued to the bottom substrate of the cell, (b) 
Images of the particle recorded at different intensity threshold levels. Red area consists of pixels that are 
used to calculate the particle’s center coordinates, (c) Apparent MSD of a glued particle at different levels 
of threshold intensity, (d) The optimum value of intensity threshold is determined from the minimum 
of MSD for the glued particle. 


so that measurement errors result in a constant additive background in MSD. The assumption of uncor¬ 
related errors almost certainly holds true for static errors, but it may not be the case for dynamic errors. 
For instance, the errors due to birefringence fluctuations are likely correlated to some extent over the 
time interval corresponding to the time scales of the fluctuations. This will add a time lag dependence to 
the background so that, in general, it is a function of t, i.e., A^. The presence of the background should be 
taken into account in data analyses. We will experimentally estimate these errors and their time depen¬ 
dence by determining the apparent MSD of immobilized particles. Digital images of colloidal particles, 
captured at a maximum frame rate of 2400 fps (time resolution is 0.4 ms), were analyzed to find the coor¬ 
dinates (x, y) of the particle’s center using the intensity-weighted algorithms. Specifically, we computed 


X = 


LiljXi 

Lih 


£ i h y i 
Liii ’ 


(4) 


where x/ and y; are coordinates of i-th pixel, 7/ is its intensity. On each frame with 8-bit gray scale (0- 
2SS arbitrary units of intensity) we take into account only the pixels that have an intensity higher than 
certain threshold intensity, 7,_th- The different threshold intensity results in different apparent MSD for 
an immobilized particle. Images of the particle at different threshold levels are demonstrated in figure]^ 
Red area consists of pixels that are taken for the calculation of the particle’s center coordinates. The 
optimal level of the threshold which gives minimum value of the apparent MSD of glued particle was 
extracted from the dependence of MSD vs intensity threshold, figurej^ 
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To establish the limit of accuracy iit the measurements of particle’s positions that depends on bire¬ 
fringence, we used particles immobilized (by a Norland adhesive) at the bottom plate of the cell hlled 
with three different fluids: two types of a nematic and water as isotropic fluid. The probing light beam 
traveled twice through the entire thickness of the ceh. The apparent mean square displacement of the im¬ 
mobilized particles vs time lag is shown in figurej^ The apparent displacements represent a cumulative 
effect of errors in measuring the particle’s position caused by the optical system of the microscope, vi¬ 
brations and birefringence. It grows with birefringence of the material, being the largest for the nematic 
pentylcyanobiphenyl (5CB) with the highest birefringence (~ 0.2). In all cases, the apparent MSB in the 
time range of interest was about 10“^® m^ or less; these values are about 100 times smaher than the MSB 
of free spheres experiencing on these timescales an anomalous diffusion described in the main text. 



Figure 5. (Color online) MSB versus time lag for immobilized silica spheres of diameter 5 pm in four 
different cells (aU of thickness 50 pm) filled with two types of a NLCs and water; the data for an empty 
cell are labeled as “dry particle”. Normal surface anchoring. Higher birefringence results in a higher 
apparent displacement of the particles. 


3. Results 

The measured MSB vs. time lag t dependencies for ci = 5 pm silica spheres in IS-8200 are presented 
in figure]^ (a), for perpendicular anchoring. In the isotropic phase (elevated temperature T - OO^C) the 
diffusion is normal in the entire range of time lags, with the diffusion coefficient D-9.2- 10“^® m^/s. In 
the nematic, at T = 50°C, the diffusion becomes anisotropic, with MSB being different when measured 
parallel and perpendicular to Hq. At relatively long time scales, x > (20-40) s, both MSB components 
grow linearly with T, with the diffusion coefficients Dy = 1.9-10“^® m^/s and Dj_ = 1.4-10“^® m^/sforthe 
normally anchored spheres. At the times shorter than about (20-40) s, the MSB time dependence becomes 
markedly nonlinear, hgure[^(a). An apparent MSB of particles glued to the bottom of the cell (see hgure|^ 
is much smaller than the MSB in the nonlinear range and is practically time-independent; it merely adds 
a constant background to the MSB of the moving particles. Thus, the non-linear behavior at short times 
is not a spurious effect due to the finite accuracy of measurements l37l or birefringence. 

To obtain a better insight into the different diffusion regimes and the characteristic times limiting 
their borders, we calculated the velocity autocorrelation function Cy||(T) = (vx{t)Vx(0)) (SSlIlll. where Vx 
is the translational velocity of the particle along the x-axis, and a similar quantity Cv± (x) - (vy (x) Vy (0)) 
for y-direction. In fact, one calculates the autocorrelation function of the mean velocity over hnite time 
interval between the position of the particle (time lag) which is much shorter than the correlation time. 
Under this condition, the mean velocity autocorrelation serves as a good estimate of the velocity autocor¬ 
relation function l40l . 

For the diffusion in the isotropic fluid and for the diffusion at large time scales in the nematic, Cv\\ (x) 
and CvL (t) are close to zero, figure (b), as it should be for the normal diffusion with MSB growing 
linearly with x. However, both Cy\\ (x) and Cpy (x) become negative in the nematic, when the time lag 
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Figure 6. (Color online) (a) MSD versus time lag for 5 ^im silica particle with normal surface anchoring 
diffusing in the isotropic (T = 60°C) and nematic (T = 50°C) phases of IS-8200, in the directions parallel (x) 
and perpendicular (y) to the overall director. CeU thickness 50 pm. (b) Velocity autocorrelation function 
calculated from experimental data for the normally anchored sphere moving in the nematic, parallel and 
perpendicular to the director and in isotropic phase. Arrows show Tsup and TjuJ), where Cp (t) is close 
to 0. 


is within a certain interval Tsup < t < Tsubi indicating o: < 1, i.e., subdiffusion. At shorter time scales 
T < Tsup. Cy|| (t) and Cp±_ (t) become positive, indicating a > 1 and superdiffusion. Determination of Tsup 
is straightforward, as a point where the velocity autocorrelation functions change their sign. The value 
Tsub is determined approximately when the deviation from zero exceeds 10“^^ m^/s^, a typical scatter of 
VACF data in figure[^(b) as it was proposed in (26l. 


4. Discussion 

The anomalous Brownian motion is only observed in the nematic phase. As soon as the nematic is 
brought into isotropic state, the sub- and superdiffusive behavior changes to a normal linear MSD time 
dependence down to the shortest experimentally accessible times. Thus, anomalous dynamics must be 
related to the dynamics of additional degrees of freedom that exist in the nematic, but not in the isotropic 
phase, namely to the nematic director dynamics. The experiments demonstrate that the orientationally 
ordered environment influences the Brownian motion of a particle most profoundly, causing, in addition 
to anisotropy, anomalous super- and subdiffusion. The corresponding times, Tsub and Tsup, vary with the 
type of anchoring at the particle’s surface, size, and displacement direction (26). Above these time scales, 
the diffusion becomes normal (but still anisotropic). The anomalous character of the diffusion in the 
range t < Tsub is evident not only in the nonlinear dependency of MSD on t, figure (a), but also in the 
behavior of velocity correlation functions, figure[^(b). 

The current models of Brownian motion in a nematic (l6l|20l[^[23 consider the director field around 
the particle as being stationary. The predicted diffusion is always normal albeit anisotropic. Our results 
agree well with these models if t is large, t > Tsub- At t < Tsub, however, the diffusion becomes anomalous; 
we attribute the effect to the director fluctuations. 

The relaxation times of director fluctuations relevant to the Brownian motion are expected to be 
limited by the interval t^ < t < where t^ = rjcf-ln^K and t/, = rjh^ln^K are the characteristic times 
associated with the particle diameter and the cell thickness, respectively. At times t > t/,, the fluctuations 
are suppressed by the surface anchoring at the cell plates. For t « t^, the influence of perturbations with 
the wavelength much shorter than d averages to zero over the distance d. For d-S pm in an IS-8200 cell 
of thickness 50 pm, the limits are estimated as =; 0.3 s and T/, 30 s. The time range < t < thus 
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embraces the experimeittally determined Xjub = (20 - 42) s and Xjup = (4 - 10) s, see figure]^ 

The nematic is a viscoelastic medium, in which the director field n(r, f) is coupled to the the velocity 
field v(r, t). Bofh n(r, t) and v(r, t) are perturbed by the particle and by the director fluctuations. Trans¬ 
lational motion is coupled to the orientational dynamics of n(r, t). In its turn, director reorientations 
induce torques and forces that cause the material to flow (the so-caUed backflow effect (^|4T]|42l) and 
thus modify v(r, t). The director fluctuations establish an intrinsic memory at the scales < r < t/,, 
which is typically much longer than the hydrodynamic memory time of the isotropic fluid (33 . The in¬ 
trinsic memory is known to cause both superdiffusion and subdiffussion, sometimes just by varying the 
parameters of the very same system @3- Below we present qualitative effects that help to understand 
the connection of director fluctuations to the intrinsic viscoelastic memory and anomalous diffusion in 
the interval t^<t < 

The equation of motion for the director fluctuations 5n in a bulk nematic liquid crystal reduces to the 
torque equation for viscous and elastic torques (33 , 

d o 

j^—6n^KV^Sn, (5) 

where yi is the rotational viscosity and K is the effective elastic constant. Performing Fourier-Laplace 
transform of equation (3 yields the dispersion relation with purely imaginary frequency w = 
from which it follows that the fluctuation modes are overdamped and thus purely relaxational with the 
relaxation time Xq = jiUKq^]. Their power spectrum is then 


Iq (ej) = 


k^T Xq 

^Kq^ l + [cjT^f 


( 6 ) 


with the corresponding correlation function 


Cq(x) = (dn_q (0)5nq (x)) = 



(7) 


Fluctuation amplitude is thus ^|5nq|^^ = Cq(0) = k^TI [Kq^). 

If a nematic is confined in a flaf cell, then, due to restrictions imposed by the boundary conditions 
and particle’s size, only a discrete set of fluctuation wave vectors is allowed. Consequently, the relaxation 
times of different fluctuation modes depend on the size of the inclusion, cell thickness and anchoring 
conditions at the surfaces I36II44I . 

A coupling between the LC director and particle dynamics may result in a variety of scenarios of a 
particle movement. A dipolar inclusion behaves as an elastic dipole with the dipole moment P - aR^ 
(fl = 2.04, reference (30)) and, therefore, interacts with inhomogeneities of the director field that arise 
due to thermal fluctuations. The dipole energy is then U - -AnKPSJn, so that the particle experiences a 
force F = -VfJ = 47iKPW (Vn). Neglecting the inertial effects, the particle then moves with a velocity v 
that is proportional to the force. 


F 

GnqR 


2KP 

^qR 


V(V-dn). 


( 8 ) 


Following the Fourier-transformation into reciprocal space, the director fluctuation component driq with 
wave vector q results in a particle velocity component 


2KP 


(9) 


Particle velocity autocorrelation function is then 


I 2KP\^ 

Cv,q(x) = (v_q(0)-Vq(x)) = ^-g^j q^(q-dn_q(0) q-dnq(x)). 


( 10 ) 


For the fluctuation mode that involves splay deformation, 5n is in the (n, q) plane and perpendicular 
to n (32), so that q-dn = qSnsinO, where 6 is the angle between q and n. Particle velocity autocorrelation 
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function is then proportional to the director autocorrelation function of equation (7), 

2 

q^sin^6(dn_q(0)dnq(T)) = (11) 

where we substitutesd the director correlation function from equation (t) and absorbed all constant pref¬ 
actors into A - {2PI3T]R]^k^TK. The obtained equation describes the contribution of thermal director 
fluctuations with wavevector q to the particle velocity autocorrelation function. To obtain a full velocity 
correlation function, the equation should be integrated over q. (Note that fluctuations with different q are 
uncorrelated.) Integration should only involve fluctuations occurring on length scales large compared to 
the particle size d - 2R, which corresponds to the wavenumbers smaller than q^-CId, where C is a 
constant of the order one. (The high-q fluctuations with the correlation length smaller than the size of 
the particle will exert uncorrelated forces on different parts of the particle, averaging out to zero.) Thus, 

C-vM-A J i7^sin^0e“^^’^‘idq, (12) 

where dq is the volume element in q-space. 

The MSD is expressed through the velocity autocorrelation function as follows I39II451 : 

T t' 

(Ax2(t)> = 2 J dt'j Cyjt]dt. (13) 

0 0 

Conversely, Cy^(T) = 1^2^ (Ax^(t)). The couphng mechanism discussed above evidently leads to 
Cy(T) > 0, and, therefore, to superdiffusion. Thus, direct dipole coupling to the director fluctuation dy¬ 
namics may be responsible for the superdiffusion that we observe. 

If the director field n(r) around the particle is bent by an external action or thermal agitation, it 
will exert a torque on the particle, proportional to the director rate of change h(t). In response, the par¬ 
ticle will rotate with angular velocity co oc h(f). Neglecting the inertial effects, a particle rotating with 
an angular velocity oj will also move with a velocity u proportional to \ld\, and thus to |h(f)|, see, e.g., 
references Thus, the particle velocity is coupled to the director fluctuations of the surround¬ 

ing nematic, and thus the particle velocity autocorrelation function is proportional to the (q-dependent) 
director angular velocity correlation function, 

Cv,q(T) = c((5h_q(0) (5hq(T)) , (14) 

where c is a coupling constat. One needs to integrate over q, noting as before that only fluctuations with 
q< qd are relevant. The particle velocity correlation function then becomes 

CvW-c J (dh_q(0) 5hq(T))dq. (15) 

Recalling that, for any mechanical property A that is a function on the phase space of a classical many- 
partlcle system, there holds (^(0)A(f)> - (j4(0)A(f)> (13, the director angular velocity correlation 

function Cn,q(T) = (h_q(0) hq(T)) = -^Ci,,q(T), where Ci,,q(T) is given by equation (7), so that 

Cri,q(T) = (h_q(0) hq(T)) = (16) 

A q Tq 

Clearly, the director angular velocity autocorrelation function is negative. This is easy to understand in 
view of the fact that director fluctuations are only small angular excursions from the mean, so that if the 
director rotates in certain direction at a given instant of time, at a later time it should be rotating back, 
which means a negative angular velocity autocorrelation. At short times, however, the relation equation 
is not valid, as the initial value of an autocorrelation function must be positive. At short times, the 


Cv,q(T) = 


I2KP 
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inertial effects cannot be neglected, and correct asymptotic behavior of Cn.qCr) is such that its initial 
second derivative is negative, so that Ch,q(T) is positive at short times, as it should be. 

From equations (Ts) and (l^ it follows that. If the particle motion couples to the director rotations, 
then the particle velocity autocorrelation function becomes 

Cv(t) = -B J q^g-^/’^qdq, (17) 

'7<'7rf 

where B is a constant. Thus, Cv(t) due to this coupling mechanism is negative, and corresponds to sub¬ 
diffusion. 

The coupling mechanism between the director bending modes and particle diffusion, discussed above, 
is indirect in the sense that director reorientations are directly coupled to particle rotations, and the lat¬ 
ter may then couple to translations. Another mechanism that would couple the director bending modes 
with particle translation is through the backflow effects (47), whereby director fluctuations induce flows 
that affect the embedded particles. In the simplest reahzation in two dimensions, the force exerted on 
the fluid and that generates a backflow depends on the director angular velocity h. Its gradient Vh, and 
the director field gradient (48). This force generates a viscose flow with a velocity v proportional to the 
force, and thus depends on h and Vh. Assuming that an embedded particle just follows the flow, its veloc¬ 
ity autocorrelation function will thus depend on the director angular velocity autocorrelation function, 
which leads to subdiffusion, as discussed in the preceding paragraph. This is most readily seen for one of 
the contributions to the backflow, which is proportional to the projection of the director angular velocity 
gradient onto the director n, i.e., nVh (48). The corresponding contribution to the particle velocity is then 
V = fcnVh, where fc is a proportionality factor. Transforming to the reciprocal space, 

iZq = fc(n-iq) hq = ifcqcosO hq, (18) 

where q is the wavevector and 6 is the angle between q and n. Velocity autocorrelation function is then 

Cy,q(T) = (u-q(O)iyq(T)) - (h_q(0)hq(T)) . (19) 

Substituting (h_q(0)hq(T)) from equation (T^ and integrating over q to obtain Cy(T), one ends up with 
the same integral as in equation (of course with a different pre-factor). Thus, the coupling to the di¬ 
rector fluctuations through backflow effects leads to a negative particle velocity autocorrelation function 
of the form of equation (17), and thus to subdiffusion. 

The described mechanisms and their interplay may lead to intricate scenarios of coupled particle dy¬ 
namics, depending on the time scales and relative strength of the two coupling effects. The time scales 
may in fact be well separated. Indeed, there are two Independent director fluctuation modes that in¬ 
volve different types of director distortions (32). If the elastic constants pertaining to these deformations 
are much different, then these two modes will relax on different time scales. In particular, only one of 
the fluctuation modes involves a splay deformation, which is the only distortion that directly couples to 
elastic dipoles through the first mechanism discussed above. 

The diffusion dynamics shown in figure [^is a complex and multidimensional process that involves 
different factors and mechanisms. This complexity does not allow one to take into account all aspects 
of the diffusion process; however, our analysis gives a physical picture where the elastic interaction be¬ 
tween the embedded particles and director fluctuations may describe the complex behavior of particle’s 
dynamics in a NLC. 

Figure]^ shows the experimental data of IS-8200 for MSD parallel to the nematic director, together 
with theoretical curves. The data are the same as in figure[^(a), with the background due to position deter¬ 
mination errors (figure subtracted. In order to generate the curves In the figure, we used the viscosity 
77 = 2.5 Pa-s (manufacturer supplied data) and the elastic constants ATu = K 33 - 30 pN and J<C 22 = 10 pN 
(our estimates), from which the average elastic constants are kTi = ^^(33 - 1 - liTn = 30 pN for splay (super- 
diffusion) and K 2 - ^K 33 + \K 22 ~ 17 pN for bend (subdiffusion) (32). Contributions of the super- and 
subdiffusion to MSD were computed by integration, equation (T^, of the velocity correlation functions 
of equation (T^ and (17), respectively. We assumed that - C/d, where C is a constant of the order 
of one. To obtain the fit in the figure, we varied C and the magnitudes of the super- and subdiffusive 
contributions. Fit results are seen to agree weU with the experimental data, both In the regions of super- 
(inset in figure[7) and subdiffusion (figure[7). 
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Figure 7. (Color online) MSD parallel to the nematic director for 5 pm particle in lS-8200 along with 
theoretical curves for anomalous (green dash line) and normal contributions (black dash line) as well as 
fit (blue solid line) of the experimental data obtained as a combination of the previous two. Inset: data in 
double logarithmic scale in order to show fittings for short time scales. 


5. Conclusions 

We investigate the thermal motion of colloidal particles in a nematic liquid crystal for the time scales 
shorter than the expected time of director fluctuations. At long times, compared to the characteristic time 
of the nematic director fluctuations, we observe a typical anisotropic Brownian motion with the MSD lin¬ 
ear in time and inversely proportional to the effective viscosity of the nematic medium. At shorter times, 
however, the dynamics is markedly nonhnear and exhibits subdiffusive (slower than oc r) evolution of 
MSD. We present a simple illustration of how the director fluctuations influence the Brownian motion 
through the long-range interactions. Although this study dealt with standard liquid crystals, the observed 
anomalous diffusion is expected to arise in any dispersive environment with reduced symmetry and 
orientational order. This is a remarkable finding, since local orientational order should profoundly in¬ 
fluence the dynamics of many complex biological systems, such as cell membranes, the cycloskeleton, 
assembles of anisotropic particles, etc. A detailed theoretical description of the effect is highly desirable. 
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AHOMa/ibHi/iii 6poyHiBCbKi/iii pyx ko/ioiahoY nacTi/iHKi/i b 
H eMaTi/iMHOMy cepeAOBi/iii^i: Bn/ii/iB (|)/iyKTyai4iii Ai<^peKTopa 

T. TypiEp® 0 . Bpofli/ii-P^, B.r. HasapeMKcP 

^ iHCTwryT c|)i3UKi/i HAH VKpaTHU, npocn. HayKH, 46, 03028 Khi'b, VKpaiHa 

^ iHCTi/iTyr piflKMX KpwcTa/iiB, KeHTCbKi/iu flep>KaBHi/iM yniBepcuTeT, Byn. VHiBepci/iTeTCbKa EcnnaHafla, 1425, 
44242 Kent, Oraiio, CLilA 

^ HapiOHa/ibHufi lexHinHi/iM yHiBepcMTeTVKpaTHU "Kni", npocn. nepeMori/i, 37, 03056 Khib, VKpaiHa 

8k 6yj\o HeutOflaBHO 0ny6./iiK0BaH0 b [TurivT. et al.. Science, 2013, 342, 1351], c|);iyKTyaL(M b opieHTapM piAKO- 
KpMCTa/iiuHoro (PK) AnpexTopa Moxtyib nepenocuTi/i Iwny/ibC bIa PKao KO/ioi'Aa, toaI A^i4>y3i^ caMoro KO.fioTAa 
crae aHOwa/ibHOK) na KopoiKi/ix npoMi>KKax nacy. BnKopncTOByronu mbtoah BiASOMiKpocKoniT ra OAHOnacruH- 
KOBOro BiAcrexeHHB, mi/i AOC.fiiA>KyBa./ii/i BunaAKOBi/iM ren/iOBi/iPi pyx KO.fioTAHoT nacTi/iHKi/i b HeMari/iuHOMy PK na 
nacax, Kopormnx 3a oniKyBannPi nac c|);iyKTyaL(iti A^^pexTopa. Ha Aoernx nacax, b nopiBHani 3 xapaKiepnc™- 
4HI/IM nacoM pe/iaKcapM HewaTMUHoro Anpexiopa, mm cnocTepira.fiM iMnoBMM aniaoTponHMM 6poyHiBCbKMM pyx 
3 cepeAHiM KBaAparoM 3MiuteHHn nacTMHKM nponoppiMHMM ftp nacy t ra oSepneHO nponoppiMHMM pp ec|)e- 
KTMBHOI B'n3KOCTi HeMaTMUHOtO CepBAOBMU^a. OAHaK Ha KOpOTUJMX Hacax, AMHaMixa nOMilHO BiApiSHBCTbCn bIa 
niHiMHoT, cepeAHiM KBaApar 3MiuteHHB b pbOMy BMnaAxy apocrae noBmbHiuje (cySAMcfjyain) a6o ujBMAUje (cy- 
nepAMtjjyaifl) bIa t. Peay/ibiaiM noncHioiOTbcn b KOHieKTCi 3B'n3Ky AMHaMixM ko./ioiahoi nacTMHKM 3 AMHaMixoio 
(()/iyKTyaL(iM AMpexTopa. 

KnioHOBi cjioBa: we/warMWHMM pipKMi/i Kpi/tcraji, 6poyHiBCbK^iCi pyx, (p/iyKTyapi)'fii/ipeKTopa 
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